Heat conduction across molecular junctions between nanoparticles 
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We investigate the problem of heat conduction across a molecular junction connecting two 
nanoparticles, both in vacuum and in a liquid environment, using classical molecular dynamics 
simulations. In vacuum, the well-known result of a length independent conductance is recovered; its 
precise value, however, is found to depend sensitively on the overlap between the vibrational spec- 
trum of the junction and the density of states of the nanoparticles that act as thermal contacts. In 
a liquid environment, the conductance is constant up to a crossover length, above which a standard 
Fourier regime is recovered. 
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INTRODUCTION 

Since the seminal work of Fermi, Pasta and Ulam [lj , 
the problem of heat transfer in low dimensional systems 
has attracted much attention (sec, e.g., Rcfs. 0,0] for a 
review), in particular with regards to nanoscale systems. 
While the problem is fundamentally important and inter- 
esting, understanding heat transport in nanoscale sys- 
tems (including molecular junctions) is of considerable 
importance for applications, e.g., in nano-electronic as 
well as thermo-electric devices where, typically, a mini- 
mal thermal conductivity is required to convert heat in 
electric current; molecular junctions are also promising 
systems for use as thermal rectifiers 0, Q . 

The problem has been the object of numerous theo- 
retical studies, which may be classified in two categories: 
(i) simple, toy-model one-dimensional systems — typ- 
ically a chain — coupled to heat baths (B-Qs (*0 nu ~ 
merical simulations of realistic systems. The latter cat- 
egory includes, among others, molecular dynamics stud- 
ies of carbon nanotubes J9|, self-assembled monolayers 
[ill ] , molecular junctions (ll| and polyethylene chains as 
well [l2|, which have been motivated to a large extent 
by experimental studies of a number of materials, in- 
cluding nanowires 13 L .LjJ, carbon nanotubes 15Hl7| and 



diamond nano-rods 18 1. 



From a fundamental viewpoint, nanoscale heat trans- 
port is closely related to the study of the validity of 
Fourier's law at small spatial scales. In the context of 
a junction of length L subject to a temperature differ- 
ence AT, the heat current J is given by 
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where K is the thermal conductivity, which in general 
is a material property. In short junctions, the physical 
dimensions of the systems of interest is comparable to the 
phonon mean free path T, and one expects a transition 



from a ballistic regime, where k ~ L when L < T, to 
a diffusive regime when L is large and temperature is 
high; in this ~ L a where < a < 1, which may 

eventually lead to the divergence of the conductivity with 
system size for systems with vanishing concentrations of 
impurities. 

The scaling of the thermal conductivity with length 
is often analysed using the Landauer formula [l9[ com- 
monly used in mesoscopic physics: 
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where T(w) is the transmission coefficient of phonons, 
and Jl and fa are the distributions of phonons in the 
two heat baths ('left' and 'right') to which the system is 
connecting. These distributions depend on the densities 
of vibrational states and on the thermal (Bose Einstein 
or classical) statistics In the limit where the temperature 
difference AT between the two heat baths vanishes and 
if the two baths have identical statistics, the Landauer 
formula becomes consistent with the Fourier law: 
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with a conductivity increasing linearly with length: k ~ 
L. 

The dependence of the thermal conductivity on the 
system length has been extensively studied in nanowires 
and carbon nanotubes Q . There has been comparatively 
fewer studies in molecular junctions, largely because mea- 
suring heat currents across molecular junctions is ex- 
tremely difficult. The usual way to measure heat current 
across a ID object consists in connecting the object to 
two electrodes and measuring the heating Joule effect due 
to a difference of voltage between the electrodes. How- 
ever, this gives access to the thermal conductance and 
one therefore needs a model for the nanoscale junction 
to extract the conductivity. Also, the thermal contacts 
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of the junction may contribute to the Joule effect, thus 
making the determination of the conductance inaccurate. 
A promising but less explored method consists in using 
nanoparticles (NPs) as thermal contacts. Indeed, NPs of 
identical or different materials can be permanently linked 
with a few nanometers long molecule — e.g. an alcenes 
or complementary DNA strands. Such particle may then 
serve as thermal contacts, either to be inserted within 
a NEMS device for thermal transport, or manipulated, 
heated and thermally monitored using optical methods. 

In the present article, we propose a numerical proto- 
col [i(J EJ to measure the thermal conductivity of a 
molecular junction permanently linking two NPs, pos- 
sibly immersed in a liquid. We examine the steady state 
current/temperature difference for the junction, and we 
compare the conductivity in the linear regime for junc- 
tions in vacuum and in a liquid. We unveil the crucial 
role of the liquid solvent on the thermal conductivity of 
the junction with respect to its length. For a junction 
in vacuum, i.e., "dry", the thermal conductivity is found 
to increase linearly with the junction length L, in agree- 
ment with the Landaucr formula, eq|3] For a junction 
immersed in a liquid solvent, or "wet" , in contrast, we 
find a different behaviour: when the junction is short, 
the conductivity increases approximately linearly with L, 
as for the dry case, then saturates for longer junctions; 
the difference results from the "friction" of the solvent, 
which damps the low-frequency modes of the junction. 




FIG. 1: Snapshot of a wet junction (with N = 13 beads) 
between two NPs embedded in a Lennard Jones liquid. The 
region beyond the dashed line is thermostatted at Tq. 



COMPUTATIONAL DETAILS 



We have used molecular dynamics (MD) simulations to 
model a molecular junction between two NPs, embedded 
or not in a liquid matrix, as illustrated in Fig. [T] The 
NPs each contain 555 atoms cut from an FCC lattice, in- 
teracting via the combination of a simple Lennard- Jones 
interaction and a nonlinear spring that forms the so called 



FENE interaction |22|. The stiffness and maximum ex- 
tension of the nonlinear spring are /cfene = 30e/cr and 
! mH = 1.5<7, with a and e the diameter and energy of 
the Lennard Jones potential (see below). The molecu- 
lar junction, whose setup is displayed in Fig. [21 consists 
of N beads (N = 5 — 25); each bead interacts with its 
neighbours via a Hookcan spring of stiffness k, which we 
varied, and equilibrium length Iq = 1.12a, exactly at the 
minimum of the Lennard Jones potential. The end atoms 
of the junction are connected each to one atom from the 
two NPs so that the equilibrium length of the junction is 
Lq = (N + 1)Iq. However, in order to keep the junction 
locally perpendicular to the NP, second-neighbour inter- 
actions are added in the form of soft springs of stiffness 
k = 30e/er and equilibrium length l' a = 1.59<r. Apart 
from the elastic springs, the atoms of the junction in- 
teract with their neighbours through the bending energy 
term H = k(1 + cos#), where 9 is the angle between 



FIG. 2: Schematic representation of the junction. The dark 
circles represent the beads forming the junction, which inter- 
act via Hookean springs (black wavy lines); the light circles 
represent NP atoms and these interact through FENE springs. 
The soft springs between the NPs and the junction are indi- 
cated by dashed lines. 



two successive segments, and k is the bending constant 
which sets the rigidity of the chain; in what follows, we 
set k = lOOOe which yields a stiff molecular junction with 
a persistence length l p ~ 2000^0 much larger than the 
equilibrium length of the junction, even for the longest 
junctions considered (N = 25). Finally, each bead of the 
junction interacts with the atoms in the NPs through the 
Lennard- Jones potential Vlj(?") = 4e((cr/r) 12 — (er/r) 6 ), 
where e and a fix the units of energy and length, respec- 
tively, which we set to 1, as we also do for the mass m of 
all atoms. For wet junctions, the particles in the liquid 
interact via Lennard- Jones potentials with the same pa- 
rameters e and <7, as do the the atoms of both the junction 
and the NPs with the liquid. In what follows, all quan- 
tities are expressed in Lennard- Jones units (LJU), which 
are based on the elementary units e, a and m. 

To prepare the wet junction in an equilibrium state, 
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time (LJ) 

FIG. 3: Heat conduction across the dry, N — 5 junction with 
k — 6000 as a function of time: energy supplied to the hot 
NP (black line), energy removed from the cold NP (red line) 
and sum of the two (blue line). 

the whole system is first equilibrated at temperature 
To = 0.75e//cs under pressure To = 0.0015e/er 2 using a 
Nose-Hoover thermostat. For the dry junction, we ther- 
malize the junction and the NPs at the temperature To 
and remove the total angular momentum of the system 
to avoid the possibility of a permanent rotational motion. 
To speed up equilibration in the latter case, we have ap- 
plied a small friction force to all atoms in the system, 
running over 500,000 time steps (dt = 0.002-y/mcr 2 /e). 
After this initial equilibration period, for both wet and 
dry junctions, one of the two NPs is heated up to con- 
stant temperature T p > T using velocity rescaling while 
the temperature of the other NP is maintained at T . In 
presence of a liquid matrix, the region at distance > 10er 
from the surface of the NP, schematically represented in 
Fig. [TJ is thermostatted at To; the whole system is also 
maintained at constant pressure To- Under these condi- 
tions, the energy transferred from one NP to the other 
across the junction corresponds to the amount of energy 
necessary to maintain the cold NP at Tq. 



CONDUCTIVITY OF DRY VERSUS WET 
JUNCTIONS 

Figure [3] displays the energy transferred across a short 
(N = 5) dry junction with k = 6000 as a function of 
time. The spring stiffness is chosen here such that the 
characteristic pulsation (angular frequency) , cu = ^Jk/m : 
matches the density of states (DOS) of the NPs, as will 
be discussed below. The junction is found to conduct en- 
ergy; Fig. [3] shows that almost immediately after heating 
the junction, the energy supplied to the hot NP is exactly 
compensated by the energy removed from the cold NP, 
and the junction enters a steady state regime: the energy 
transferred from one NP to the other increases linearly 
with time, i.e., the power T flowing through the junction 




FIG. 4: Power as a function of the temperature difference 
AT for a short (N = 5) dry junction (o, black line) and wet 
junction (□, red line). The lines are guides to the eye. 

is constant. 

The corresponding results for the power as a function 
of the temperature difference AT between the hot and 
cold NP are presented in Fig. S) As expected, the power 
increases with AT. A linear behaviour, is observed for 
temperature difference smaller than 1 LJU, which corre- 
sponds in real units to hundreds of Kelvins. This linear 
behaviour allows one to define unambiguously the con- 
ductance of the junction. With these temperature differ- 
ences, the typical power is a few tenths of nW. For larger 
AT's, the power increases more slowly, a consequence of 
anharmonic contributions. Indeed, the atoms of the NP 
are connected with non linear springs which, at high tem- 
perature, may couple strongly with the springs from the 
junction, making the hole system strongly anharmonic. 
Fig. @] shows also that the energy transfer is not affected 
by the presence of a liquid surrounding the system, at 
least for short junctions; this suggests that the heat flow 
via the liquid is negligible compared to that across the 
junction. 

To verify this statement, we have run simulations with- 
out a junction, i.e. for two NPs at a surface-to-surface 
distance equal to the equilibrium length of the N = 5 
junction considered above, To = 6/o — 7er. In practice, 
this is achieved by attaching the center of mass of each 
NP to a fixed point using a stiff spring. The two NPs are 
immersed in the same liquid as the NP-junction system 
and equilibrated, after which one of the NP is heated up 
to temperature T while the other is thermostatted to To 
(and so is the liquid at distance > 10ct). In this situation, 
the liquid between the two NPs conducts energy due to 
the local temperature gradient (which is ~ 1 K/nm!). 
The conductance of the liquid is obtained by measuring 
the energy necessary to thermalize the cold NP. 

Fig. E] provides a comparison of the energy transferred 
across short wet junctions for different values of the 
spring constant k to the energy transferred across the 
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FIG. 5: Energy transferred across the wet junction compared 
to the energy transferred across the liquid, for a short (N = 5) 
junction and for a pure liquid with a comparable interparticle 
separation, respectively d = 5 (dashed lines) and d — 6cr 
(dashed-dotted lines). The temperature difference between 
the two contacts of the junction is AT = 1.15. 



liquid alone. The range of values of k has been chosen so 
as to cover all possible situations of overlap between the 
DOS of the NPs and the vibration spectrum of the junc- 
tion. In particular, for the lowest value (k = 100), the 
characteristic pulsation uj = yjk/m is smaller than the 
low frequency branch of the NP DOS, and one therefore 
expects weak energy transfer in this case. Indeed, the 
energy transferred across the liquid is negligibly small 
compared to the energy transferred across the junction. 
The power flowing across the 1 nm thick liquid layer is 
approximately 1 nW. 



So far, we have considered a perfectly wetting (hy- 
drophilic) situation. For large contact angles (i.e., a more 
hydrophobic particle) — which can be achieved by tuning 
the attractive part of the NP-liquid interaction — one ex- 
pects even smaller energy transfer capability across the 
liquid. Indeed, at the interface between the NPs and 
the liquid, the temperature field displays a discontinuity 
due to the finite conductance of the interface [2^]; this 
discontinuity is known to increase with hydrophobicity, 
as the mismatch between the two media decreases the 
transmission of phonons through the interface. For our 
two NPs and given the length L of the junction, the ef- 
fective temperature gradient in the liquid decreases as 
hydrophobicity increases because the temperature jump 
at the NP-liquid gets larger (i.e., the temperature of the 
liquid is smaller). Hence, the effective conductance of the 
liquid is expected to be even lower than in the hydrophilic 
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FIG. 6: Power as a function of spring stiffness k, for wet (□) 
and dry (o) short (N = 5) junctions; AT = 1.15. 



INFLUENCE OF THE STIFFNESS 

We now examine the variation of the conductance with 
the stiffness of the springs connecting the atoms within 
the junction. The results are presented in Fig.[6]for both 
the wet and the dry A" = 5 junction. The dependence on 
k is found to be non monotonous in both cases: the power 
first increases, reaches a maximum — at the same value 
of k (~ 2000) for both wet and dry junctions — , then 
decreases, the dry junction thus appearing to be a better 
conductor than the wet junction for a large stiffness. 

This particular behaviour can be explained qualita- 
tively in terms of the DOS of the NPs and the vibrational 
spectrum of the junction, as depicted in Fig. [7] The DOS 
of the NPs has been obtained by Fourier transforming 
the velocity autocorrelation function. For the junction, 
we have only considered, for the sake of clarity, the spec- 
trum of stretching modes, whose resonance frequencies 
are given by uj p = y/k/m sin(7rj> /2-/V) , p being an integer 
between 1 and N— 1. Apart for these A^ — 1 modes, there 
arc N — 2 resonance frequencies corresponding to bend- 
ing which, in the harmonic approximation, are decoupled 
from longitudinal modes so that the corresponding spec- 
trum does not depend on k. Different cases may be en- 
countered, depending on k, as demonstrated in Fig. m 
For small values of k, only a few modes of the junction 
overlap with the low frequency branch of the NP DOS, so 
that only the low frequency modes of the NP are excited 
and a small conductance is expected, as indeed observed 
in Fig. |U When k ~ 2000, the matching between the 
stretching modes and the NP DOS is optimal, and the 
junction is a very good conductor, consistent with the 
results of Fig. [S] For large values of k, the normal modes 
of the junction have high frequencies and there is little 
overlap with the NP DOS so that the energy transfer 
across the junction is expected to decline. 

It turns out, however, that the conductance of stiff 
junctions is larger than that of soft junctions, cither dry 
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FIG. 7: Vibrational spectrum of a TV = 5 junction (vertical 
lines) superimposed on the NP DOS (full line) for different 
values of the longitudinal stiffness k: k — 200 (red lines), 
2000 (black lines) and 60000 (dashed line in this case, only 
the lowest frequency mode is shown). 




FIG. 8: Conductivity of dry (o) and wet (Exjunctions as a 
function of length, measured by the number of beads of the 
junction N. Here, k = 6000 and AT = 1.15. 

or wet (cf. Fig. [6J . This may appear counterintuitive, as 
there is some overlap between the soft -junction spectrum 
and the NP DOS while there is none for the stiff junc- 
tion. A possible explanation is that, for large stiffness, 
the high-frequency longitudinal modes of the junction ex- 
cite higher harmonics of the NPs, allowing energy to be 
transferred between them. As noted earlier, Fig. [6] also 
reveals that the dry junction conducts heat better than 
the wet junction. The difference may be understood in 
terms of the friction of the liquid, which damps the low 
frequency modes of the NPs, thus reducing the efficiency 
of the non-linear coupling discussed above. 



EFFECT OF LENGTH 

We discuss, finally, how the length of the junction af- 
fects the energy transfer between the NPs. The variation 



of the conductivity of dry and wet junctions with length 
is displayed in Fig. [H The stiffness k (= 6000) has been 
chosen here so as to ensure a good matching between 
the junction and the thermal contacts with the NPs. 
For the dry junction, the relation between conductiv- 
ity and length is monotonous: the conductivity increases 
almost linearly, i.e., the conductance is approximately 
constant. This behaviour is standard for unidimensional 
systems, and easy to understand. Upon increasing the 
length of the junction, the number of vibrational modes 
iV mot jes = 2iV — 3 increases proportionally and the overlap 
with the NP DOS increases, provided the normal mode 
frequencies w str otching = \/ k/m and ^bending = V 'k/1% 
are not too small. This is demonstrated in Fig. [S] where 
we superimpose the spectra of free dry junctions (not con- 
nected to nanoparticles) of two different lengths (N = 5 
and N = 20) with the NP DOS: the overlap clearly in- 
creases with length and, as a result, the conductivity in- 
creases. 

For wet junctions, the scenario is different, as can be 
seen in Fig. HI The conductivity first increases with 
length, following closely the results for the dry junction; 
for N > 13, however, the current starts decreasing with 
increasing length, i.e. the conductance drops. This non 
monotonous behaviour may be understood by analysing 
the vibrational spectra of wet junctions, also shown in 
Fig. IH1 These have been calculated from simulations of 
a free junction (i.e., without NPs) immersed in a liquid. 
For the shortest junction (JV = 5), the DOS of wet and 
dry junctions are very similar, except perhaps at low fre- 
quencies; this is evidently consistent with the fact that 
the conductivities of short junction is not affected by the 
presence of the surrounding solvent, as observed above 
(Fig. [51). Fig. [9] however shows that the spectra of long 
wet and dry junctions may be different and, in particular, 
the low-frequency modes of the dry junction are absent 
in the wet junction. We infer from this observation that 
the surrounding solvent damps the low-frequency modes 
of the wet junction. To ascertain this, we also show in 
Fig. [9] the DOS of the Lennard- Jones liquid, here ther- 
mostatted at the reference temperature T = 0.75. It 
is found to be significant only at low frequencies — u> 
less than 40 or so; in contrast, for higher frequencies, the 
spectrum of the liquid is vanishingly small. Hence, high- 
frequency modes of the junction feel the presence of the 
solvent as a frictionless elastic medium, and are thus not 
affected by the presence of the liquid. This explains why 
the high-frequency part of the spectrum of dry and wet 
junctions are very similar, irrespective of the length of the 
junction. Low-frequency modes, on the other hand, feel 
the solvent as a viscous medium, and are therefore signif- 
icantly damped by the solvent. This is can in fact be seen 
in Fig.[9]where the low-frequency modes are almost com- 
pletely absent, more precisely almost completely damped 
— in the spectrum of the long wet junction. 
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FIG. 9: Vibrational spectra of dry (black lines) and wet (red 
lines) junctions with N = 5 (top) and TV = 20 (bottom), 
superimposed on the DOS for the NPs (blue solid lines) and 
that of the liquid (dashed blue lines). 



CONCLUSION 

We have studied the problem of heat transfer accross 
a molecular junction between two nanoparticles that lie 
either in vacuum and in a liquid environment using clas- 
sical molecular dynamics simulations. While our model 
remains "generic" , it does capture the essential physics of 
the problem. For the junction in vacuum, the well-known 
result of a conductance that is independent of length is 
recovered; however, the conductivity depends sensitively 
on the overlap between the density of states of the junc- 
tion itself and that of the nanoparticles to which it is 
connected. The presence of a liquid is found to affect 
this behaviour in a significant manner: the conductance 
is constant to a crossover length, above which it starts 
to decline and the standard Fourier regime is recovered. 
Our results can be rationalized in terms of the overlap be- 
tween the vibrational spectrum of the junction and that 
of the immersing liquid, which is qualitatively different 
in short and long junctions. 
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